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We calculate the effect that a primordial homogeneous magnetic field, So, will have on the different 
CMB power spectra due to Faraday rotation. Concentrating on the TB, EB and BB correlations, 
we forecast the ability for future CMB polarization experiments to constrain Bo- Our results depend 
on how well the foregrounds can be subtracted from the CMB maps, but we find a predicted error 
between = 4 x 10~^^ Gauss (for the QUIET experiment with foregrounds perfectly subtracted) 
and 3 x 10^^"Gauss (with the Clover experiment with no foreground subtraction). These constraints 
are two orders of magnitudes better than the present limits on Bo. 



I. INTRODUCTION 

Magnetic fields, with amplitudes of around 1/iGauss 
have been observed in galaxies and clusters [1-5]. Indeed 
it has been argued that these fields may play a role in 
galaxy formation and evolution. It is believed that these 
fields may have arisen from a common mechanism: an 
amplification of a weak primordial field through adia- 
batic collapse or some form of cosmic dynamo. There is 
a suite of proposals for the origin of the primordial field, 
i.e. the seed field, from phase transitions, to parity vio- 
lating processes and inflation. In general the seed fields 
generated by these mechanisms will be stochastic and 
spatially varying. There is, however, the possibility that 
the primordial magnetic field is homogeneous, embedded 
in the large scale fabric of the universe [6- 12] . All of these 
possibilities should lead to distinct observational signa- 
tures. In this paper we will focus on one such signature, 
the polarization of the Cosmic Microwave Background 
(CMB). 

The polarization of the CMB is generated by Thomson 
scattering of the quadrupole temperature anisotropy at 
the last scattering surface. Since the dipole anisotropy 
does not generate any polarization, the amplitude of the 
polarization signal will be significantly lower than for the 
temperature anisotropics, which get most of the power 
from the scattering of the dipole anisotropy. This low 
amplitude leads to experimental challenges. However, 
the different polarization power spectra will provide im- 
portant and complimentary information to what we can 
learn from studying the temperature fluctuations alone, 
and up and coming experiments, such as Clover [13, 14], 
EBEX [15, 16], QUIET [17, 18], Spider [19] and Planck 
[20, 21], will attempt to do measure the CMB polariza- 
tion signal with unprecedented precision. 

Primordial magnetic fields leads to Faraday rotation of 
CMB polarization. This alters the overall pattern of the 
CMB polarization by mixing E and B modes. The effect 
can play a role on small scales (in the case of stochas- 
tic magnetic fields) [22-26] and large scales (in the case 
of homogeneous magnetic fields) [27-30]. In this paper 



we explore the latter and forecast how well we can use 
future CMB polarization experiments to constrain the 
amplitude of the magnetic field. 



II. CMB POWER SPECTRA IN THE 
PRESENCE OF A HOMOGENEOUS MAGNETIC 
FIELD 

The polarization of the CMB can be described in terms 
of the Stokes parameters, Q and U . These parameters 
quantify the intensity of the incoming radiation in mu- 
tually orthogonal directions. The Q and U parameters 
are not rotationally invariant and depend on the specific 
choice of coordinate system. It is therefore convenient to 
transform to rotationally invariant E (curl free) and B 
(divergence free) modes before comparing observations 
to theory. These are defined by 



{Q±iU){n) 



{Elm ± iBlm)±2Ylmin), (1) 



where ±2^m('^) denotes the spin- weighted spherical har- 
monics. 

The Thomson scattering process will only generate E 
modes from scalar perturbations of the metric and stress 
energy tensor. In the standard cosmological picture to- 
day the B modes are only generated by relic gravita- 
tional waves from infiation or from gravitational lensing 
of E modes. These two processes will induce a frequency 
independent B mode signal for low Is and high Is re- 
spectively. B modes created by gravitational waves and 
lensing will lead to a non-zero BB autocorrelation signal. 
Non-zero correlations between the B mode and the E or 
T (temperature) modes can only occur under processes 
that break parity invariance. Thus, any detection of a 
non-zero TB or EB signal would be a sign of physics 
beyond the cosmological standard paradigm. 

A large-scale homogeneous magnetic field is an exam- 
ple of a parity violating field, and the Faraday rotation 
of CMB polarization induced by such a field will indeed 
result in non-zero EB and TB correlations. In an other- 
wise standard cosmological scenario such a magnetic field 



2 



will tend to rotate E modes generated by Thomson scat- 
tering into B modes. Thus one would expect non-zero 
BB, TB and EB signals with power spectra similar to 
the corresponding E mode spectra, but with a different 
amplitude. Since Faraday rotation is frequency depen- 
dent, one would expect the magnitude of this effect to 
depend on the frequency in which the CMB polarization 
is observed. 

The effects of such a large-scale magnetic field on the 
CMB power spectra were studied [29] and [30]. In [29] 
the authors looked at the TB power spectrum and found 
that one could expect to put upper limits on a primor- 
dial magnetic field of the order Bq ^ lO^^Gauss by the 
Planck experiment, where Bq denotes the strength of the 
magnetic field today. In [30] the authors derived equa- 
tions for all the B power spectra using the total angular 
momentum method [31]. Our method for calculating the 
different CMB power spectra is based on ref. [30], and 
we refer to this paper for further details on the notation 
and derivations. 

The effect of the Faraday rotation is conveniently de- 
scribed by the time independent parameter F, given by 



F = 0.7 



_e \ / lOGHz V 



10-9Gauss J \ I' J 



(2) 



where B denotes the strength of the magnetic field and 
z/ is the frequency of the radiation. In [30] the authors 
found that, if there are no B modes in the absence of a 
magnetic field, the effect on the CMB E and B modes 
by including a non-zero F can be expressed by 



Elm ± iBlm = ^ 



(-1)' -' + 1 



± 



(-1) 



R{lm, l',l"m)Ep^.{3) 



where R{lm,l',l"m) can be expressed in terms of 
Clebsch-Gordan coefficients: 



mm,l',l"m) = ^'^J^^ 

X {l',l";m-m",m"\l'l";l,m) 

X {l',l";0,2\l'l";l,2). (4) 

The effect of the magnetic field is hidden in the quantity 
Ep- 



Im 



X / d-qfe-^21' + l)ji,{FT)P^"\ri) 
Jo 

X el°)[fc(%-^)]. (5) 

This quantity only differs from the standard expression 
by the extra factor {21' + l)jV(Fr), where ji denotes 



spherical Bcssel fimctions of order I and t is the opti- 
cal depth. Note that limx^oji{x) = for I > such that 
all the introduced corrections reduce to for F — 0. For 
exact definitions of the other quantities appearing in (5) 
we refer to [30]. 

With this definition of Ej^^' the TE power spectrum is 
now given by 
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(6) 



Resulting TE power spectra for different values of F are 

shown in Figure 1. 

For the other power spectra we will only study the 
limit <C 1, which is likely to be a good approximation. 
In this limit we can neglect contributions from T > 1 in 
(5). The TB correlation can then be approximated by 
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Here it should be stressed that the TB correlations are 
non-diagonal. In the limit of F <C 1 the contributions 
come from the Ti±iBi correlations. The Cf^^ given in 
(7) then represents a mean of the two off-diagonal contri- 
butions. Similarly, the EB correlation can be expressed 
as 



EB 



21 + 



' r(;2_^2)(;2_4)lV2 
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(8) 

Also here the correlation is off-diagonal and eq. (8) is the 

mean of the Ei±iBi contributions. 

The autocorrelation polarization power spectra are 
given by 
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The resulting power spectra have been calculated with 

a modified version of CAMB [32, 33]. In Figure 1 the 

power spectra are plotted for different values of F with 

the other cosmological parameters fixed to concordance 

ACDM values. From the figure it is obvious that both the 

TB, EB and BB power spectra are promising targets for 

constraining Bq, especially since all these power spectra 

arc expected to be zero in the case of Bq = 0. Note 

that in the limit F < 1 Cf^ and Cf scale hnearly 

with F, such that the quantity Cf^/F is constant. The 

amplitude of the C^^ power spectrum will however scale 
likeF^ [30]. 



III. FORECASTS OF EXPERIMENTAL 
CONSTRAINTS 

The aim of this work is to forecast the sensitivity of 
future CMB polarization experiments to the effect of a 
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FIG. 1: CMB power spectra for models with F=0, 0.1, 0.5. 
The plotted quantity is 1(1 + 1)C^^ /2n as a function of mul- 
tipole I. See text for definitions of the different C^^ . 



homogeneous magnetic field as described in Section II. 
We will do this using a Fisher matrix approach. In our 
analysis we will assume all cosmological parameters ex- 
cept for F to be known. This can be justified since our 
results are driven by the TB and EB correlations, which 
are zero when F = 0. We will also make use of the BB 
correlation, but we note that the BB power spectrum 
from a homogeneous magnetic field differs significantly 
from the BB correlations from gravitational waves and 
lensing of E modes. The amplitude and shape of the 
different B power spectra will indeed depend on the pa- 
rameters governing the E modes, but these parameters 
are already known with relatively good accuracy, and the 
precision will of course improve further with upcoming 
experiments. 



A. The Fisher matrix approach 

Fisher matrices are a commonly used tool for forecast- 
ing experimental sensitivity. Here one uses the curvature 
of the likelihood function around some fiducial model to 
predict the error of theoretical parameters given informa- 
tion on experimental parameters of the experiment being 
considered. In our case, since F is the only parameter 
that is allowed to vary, the Fisher matrix will be 1 x 1, 
and the error in F is given by 



(aery 



(9) 



where X and Y are {T, E^ B} and af^^ denotes the ex- 
pected error in in a given experiment. These errors 
are given by [34, 35] 



/sky(2^ + l) 

;cr + /scan-^Vsr')' 

{Cf"" + /scan^^x'x^r') 
ipj^ -I- /scanWyyS; 
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(10) 



Here /g]jy denotes the fraction of the sky covered by the 
experiment. This factor accounts for the smaller num- 
ber of independent samples achieved when covering only 
parts of the sky. For a full-sky experiment like Planck, 
this factor is the part of the sky actually used in the anal- 
ysis after the galaxy cut has been made. The fraction of 
the sky that is scanned by the experiment is denoted 
/scan- It appears multiplying w^^ because a smaller 
scanned sky fraction means longer integration time per 
pixel. Partial-sky experiments will of course aim to avoid 
the troublesome galaxy plane in their scanning strategy, 
making /sky ~ /scan for such experiments. 

A convenient notation to describe the noise in a 
pixel independent way is 

(4^s2)/(/ 



"^^ '^Vx^pix^o ) [35, 36]. Here ?L 



f^pix^Vi 



pix 



pix 
is the angu- 



lar pixel size, is the time used observing each pixel 
and A^pix is the total number of pixels. The detector 
sensitivity is denoted by s and To = 2.73K is the average 
sky temperature. We have wxy = for X 7^ F [34]. For 
X = y we have 



= 2.14 X \Q~^H 



-1 ( 



V200AiKVsec 



(11) 



where tyr is the observing time in years. From now on we 
denote wt = wtt and wp = wee — wbb- For bolomet- 
ric experiments the polarization detectors will observe 
only the Q ov U Stokes parameter at a time, leading to 
uip^ — 2wj.^. 

In eq. (10) Bi is the assumed Gaussian experimental 
beam, Bi — e~' '^e/^, where ag is the beam width, ag — 
FWHM(beam)/\/81n2. The C/s appearing in (10) are 
given in dimensionless AT/T units. When estimating 
parameter errors from the BB power spectrum, we use 
instead of ap'^ in eqn. (9), since the BB power 
spectrum scales as 



B. Foregrounds 

In the procedure described above, we have neglected all 
effects from foregrounds. How much the foregrounds will 
affect the obtained limits on the magnetic field strength 
will depend on the specific experiments used, our un- 
derstanding of the statistical properties of the different 
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foregrounds, and whether a significant proportion of the 
foregrounds has been successfully subtracted from the 
CMB maps before the parameter analysis is performed. 
To do a proper analysis of all the involved foreground ef- 
fects is outside the scope of this work. However, we will 
use a simple foreground model to estimate the possible 
effects of foregrounds on our results. 

At the frequencies relevant here {ly < lOOGHz), the 
synchrotron emission from our galaxy is the dominant 
source of polarized foregrounds [37]. For frequencies of 
order lOOGHz the polarized emission from vibrating dust 
will also contribute. In the following analysis we will only 
include these two foreground components. We will then 
simply add the foreground signal to the CMB signal in 
eq. (10) such that 



C, 



XY 



riXY 



riXY 

'-'1 



l.dust 5 



(12) 



where G^^^mb is the CMB signal and Cf,]^„^ and Cf^^.i 
come from the synchrotron and dust foregrounds, respec- 
tively. For both the foreground components we will as- 
sume a simple power law dependence [37]: 



C 



;,{sync,dust} ~ ^i^-^> ' 



(13) 



where A is the amplitude of the unpolarized radiation and 
■p is the fractional amplitude of the polarized components 
and /3 is a constant describing the scale-dependence. 
The numerical values of A, /3 and p for our dust and 
synchrotron foreground models are taken from [37] and 
quoted in Table I. 
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Synchrotron p 
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A 




101/iK 
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1.95 1.3 


Dust p 


1 


0.0098 0.0024 


A 




24AtK 



TABLE I: Numerical values for the foreground parameters 
in equation (13) for synchrotron and dust emission. Here 
P={E,B}. 

The foreground power spectra will also be frequency 
dependent. For the synchrotron component, the fre- 
quency dependence is given by Osync('^) c< c{v)v~'^-^ , 
where c{v) = (2 sinh(a;/2)/2)2 and x = i//56.8GHz. 
Osync(i^) is normalized to be 1 for ^ = 19GHz. 

For the dust emission the frequency dependence is 
modelled by [37] 



©dust = c(t/)c,(t^) ^^^^^^^^^^ ^ 



1 



(14) 



where h is the Planck constant, k is the Boltzmann con- 
stant, and Ct,{v) cx l/x^. The normalization is such that 
©dust = 1 for 1/ = 90GHz. 



We are assuming that no foregrounds have been sub- 
tracted from the maps before the analysis is performed. 
This is a very pessimistic assumption, since the multi- 
frequency design of the different experiments will be used 
to subtract different foreground sources. Also, the up- 
coming C-BASS experiment [38] will make a 5GHz full- 
sky polarization map, tailored for synchrotron subtrac- 
tion from CMB maps. The QUIJOTE experiment [39] 
will make measurements of microwave polarization on 
several low-frequency channels to further improve our 
knowledge on the polarized galactic foregrounds. We 
therefore believe that the actual ability of future CMB 
experiments to constrain primordial magnetic fields will 
be somewhere between the obtained results with and 
without the added foreground model. 



C. Dependence on experimental parameters 

We will now turn our attention towards the different 
experimental parameters, and study how they affect the 
experimental sensitivity to Bq. We compare the limits 
obtainable from the TB, EB and BB power spectra. 

The relevant experimental parameters that appear in 
af^ are /gj^y, /scan, '"'p^ and ae- When looking at 
a certain subset of the experimental parameters, we will 
fix the others to typical values for the 30GHz channel of 
the Planck experiment (see Table II). 
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FIG. 2: The contribution to l/a'p from each multipole I rela- 
tive to the sum over all Is. The experimental parameters are 
set to Planck values. The red solid lines show the results for 
the TB power spectrum, black dashed lines for the EB, and 
blue dotted lines for BB. The narrow lines show the results 
when no foreground is added, while the broad lines are the re- 
sults when using the synchrotron and dust foreground model 
described in the text. 



The effect of changing /gj^y depends partly on which 
multipoles that contribute to the sum in eq. (9). In Fig- 
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FIG. 3: aBo as a function of /gj^y Here the rest of the 
experimental parameters are set to typical Planck values. 
For each of the power spectra, (tsq is normalized to 1 for 
■^sky ~ ./scan = 1- The labeling of the graphs is the same as 
in Figure 2 



ure 2 we have plotted for individual multipoles. It 
turns out that the main contribution to the sum comes 
from Is between 80 and 500. Comparing to the power 
spectra in Figure 1, this corresponds to the first peaks 
of the power spectra. Note that for the BB power spec- 
trum the distribution is more smeared out because of 
the dependence. We also see that the general pat- 
tern does not change much when adding the foreground 
model. From Figure 2 we would not expect that a large 
experiments multipole range would be important to place 
strong constraints on Bq. 

In Figure 3 we illustrate the dependence of aso on 
fg^y in a Planck- type experiment. Here we have set 
■^sky ~ /scan and only included multipoles / > 18O°/0 
in the sum, where 4> is the typical angular extension of 
an experiment with a given /g]jy We find significant 
changes in the behaviour with and without the fore- 
ground model added. For the foreground free model, 
the overall dependence on /gj^y is not very strong, as 

long as /sky > 10^"*. Also, we see that we get the best 
constraints from experiments with /g]jy ^ 10^^, which 

corresponds to including Z > 80. This is exactly what we 
would expect from the I dependence shown in Figure 2. 
However, when the foreground is added, the importance 
of having a large /gky is more pronounced. This can be 
understood by looking at equation (10). The gain of hav- 
ing a small /scan to obtain better sensitivity is less when 
the larger, foreground contaminated C/s are added to the 
noise term. 

Next we have looked at the effect of changing the 
beam width for a Planck- type experiment. The depen- 
dence of a Bo on the beam FWHM is shown in Figure 4. 



We see that the dependence is approximately linear for 
beam widths between and 100 arcminutes, and that 
changes roughly an order of magnitude in this range for 
the BB power spectrum, and slightly more when using 
the TB and EB power spectra. There is no big change 
in the importance of the beam size when adding the fore- 
ground model, but we notice that there is a slightly larger 
relative gain of having a smaller beam when including 
foregrounds. 
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FIG. 4: (JBp as a function of the beam FWHM. Here the rest of 
the experimental parameters are set to typical Planck values. 
obq is normalized to 1 for the Planck beam of 33 arcminutes. 
The labeling of the graphs is the same as in Figure 2 



Large effects are also achieved when varying Wp and 
w^^ . In Figure 5 we show contours of crgj, in the 
w^^ — Wp^ plane when using the TB power spectrum. 
Here no foregrounds are added. It is obvious that a good 
experimental sensitivity and/or a long integration time is 
crucial to place tight constraints on Bq. From the plot we 
see that the main importance is to have good polariza- 
tion sensitivity, but as already mentioned, Wp^ and w^^ 
are proportional in standard bolometer experiments. In 
Figure 6 we show the dependence of aso on Wp^ when 
using the different power spectra {w^^ = O.bwp^ in the 
TB case). We see that improving Wp^ is more important 
for EB than in the case of TB and BB. Since the BB 
power spectrum scales as F^, the potential for using it 
to constrain Bq will be less important as the sensitivity 
improves. Another important feature in Figure 6 is the 
large difference between the behavior with and without 
foregrounds included. We observe a significantly smaller 
gain in sensitivity to a^g by decreasing Wp^ when the 
foreground model is included. 

In Figure 7 we have plotted a^o a-s a function of the 
observed frequency. In the foreground free scenario, 
this will simply result in a dependence. Including 
the synchrotron and dust foregrounds makes the picture 
more interesting, as the amplitude of the synchrotron 
power spectrum has an even stronger negative correla- 
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FIG. 5; Contours of aBn 

/ 10"^° Gauss in the wp — wt plane 
when using the TB power spectrum. Here the rest of the 
experimental parameters are set to typical Planck values. The 
Planck values of wp and wt are marked with a red cross. No 
foregrounds are added here. 
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function of wZ^ for the TB, EB and BB 



6: (tsq as a iuiii^uiuii uu p 
power spectra. Here aso is normalized to 1 for wp = 13.2 x 
10~^^ (the Planck value). Other experimental parameters 
are set to Planck values. For the TB spectrum we have used 
= O.SiOp^. The labeling of the graphs is the same as in 
Figure 2 



tion with the frequency than the Faraday rotation effect. 
We see that there is a real competition between different 
effects at low frequencies. For an Planck like experiment, 
observing at frequencies lower than 20GHz will in fact be 
less sensitive to with our foreground model. 

To summarize, we need experiments characterized by 
good polarization sensitivity and preferably also a small 
beam width and low frequency to be able to place tight 
constraints on Bq. With unsubtracted foregrounds, the 
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FIG. 7: Here obq is plotted as a function of the observed 
frequency. The labels are the same as in Figure 2. The values 
are normalized to 1 for 30GHz. 



importance of a small beam persist, while the impor- 
tance of a small w~p^ decreases. Also, with foregrounds 
included, there is preference for a larger observed sky 
fraction than in the foreground free scenario. 



D. The experiments 

We now consider some specific upcoming CMB experi- 
ments and forecast their ability to constrain Sq- We will 
concentrate on the lowest frequency channels of the ex- 
periment, given that the effect of Bq carries such a strong 
frequency dependence. 

The Planck experiment [21] is a satellite based full- 
sky experiment to be launched in 2008. Some of the 
experimental characteristics are listed in Table If. In 
addition we will assume /g^y — 0.8 due to the galaxy cut 
that has to be made before the data analysis. We have 
considered a 14 months survey. For Planck we will use ^s 
between 2 and 2500 in our analysis. 

The QUIET experiment [17, 18] is a ground based ex- 
periment dedicated to CMB polarization, which will start 
taking data in 2008. Instead of using bolometers, QUIET 
uses arrays of HEMT based receivers that measure Q and 
U simultaneaously. The experiment is divided into two 
phases. The first phase (PI) will last for 2 years using a 
Im telescope. The following phase 2 (P2) will use a set of 
three 2m telescopes and one 7m telescope. Experimental 
parameters [40] are listed in Table II. For QUIET we will 
use ^s between 50 and 2500. 

The Clover experiment [13, 14] is another ground based 
CMB polarization experiment. The 96GHz telescope 
which we will refer to here, will start taking data in 
2009, and collect data with an effective integration time 
of ~ 0.8 years. For Clover we use Is between 20 and 1000. 

We will also look at two balloon-borne experiments; 
Spider [19] and EBEX [16]. Spider will cover about 60% 
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of the sky during ^ 25 flight-nights, starting from Aus- 
traha in December 2009. We apply multipoles between 
1 = 4 and I = 500 for the Spider experiment. We will as- 
sume /sky ~ ^-^ other experimental parameters we 
refer to Table II. EBEX will having a total integration 
time of about 14 days, and start operating in 2008. It 
will cover multipoles between 20 and 1500, operating on 
three frequency bands, the lowest one at 150GHz. Fur- 
ther experimental characteristics [41] are given in Table 
11. 

The ground based QUIJOTE experiment [39] is de- 
signed to determine the characteristics of the polarized 
galactic foregrounds, covering about 10% of the sky 
on the northern hemisphere. Unlike the C-BASS ex- 
periment, QUIJOTE is sensitive to frequencies as high 
as 30GHz, which should make it suitable for detecting 
the polarization signal from primordial magnetic fields. 
This is provided that it surveys patches sufficiently far 
away from the highly foreground dominated galactic 
plane. Experimental characteristics [42] for the QUI- 
JOTE 30GHz channel are given in Table II. 



Experiment 


i/(GHz) /scan 


beam(') 


-1 


-1 
in-i« 


Planck 


30 


1 


33 


6.6 


13.2 


QUIET PI 


40 


0.03 


41 


8.1 


0.98 


QUIET P2 2m 


40 


0.04 


23 


0.67 


0.062 


QUIET P2 7m 


40 


0.01 


9 


1.34 


0.12 


Clover 


97 


0.23 


8 


0.045 


0.090 


Spider 


96 


0.6 


58 


0.12 


0.24 


EBEX 


150 


0.01 


8 


0.18 


0.35 


QUIJOTE 


30 


0.1 


22 


1.3 


2.7 



TABLE II: Experimental characteristics for some CMB po- 
larization experiments. The beam size refers to the FWHM 
value in arcminutes. 



IV. RESULTS 

Throughout our analysis we have used as a fiducial 
cosmological model a standard flat ACDM model based 
on the mean parameter values for the WMAP 5 year 
results [43] and defined by {fli)h'^,flch'^,h,T,ns,As} = 
{0.023,0.11,0.72,0.09,0.96,2.3 X IQ-^}. The parame- 
ter definitions used correspond to the standard defini- 
tions in CAMB. Here fl^h^ and ilc/i^ are the physical 
densities of baryons and cold dark matter, respectively 
(we assume massless neutrinos). The Hubble parame- 
ter is parametrized by Hq = lOOft-km s~^Mpc, and r is 
the optical depth at reionization. The primordial power 
spectrum is parameterized by the amplitude Ag (at 0.05 
Mpc"'^) and the tilt Ug. We do not expect our results 
to change significantly by minor changes in this fiducial 
model. 

The resulting forecasts for for the different exper- 
iments are summarized in Table III. 



Experiment 


TB 


EB 


BB 


Planck 


2.1(2.5) 


8(12) 


15(18) 


QUIET PI 


1.1(3.6) 


0.7(3.1) 


3.3(10) 


QUIET P2 2m 


0.15(1.2) 


0.07(0.8) 


0.5(4.3) 


QUIET P2 7m 


0.10(1.0) 


0.041(0.7) 


0.31(3.7) 


Clover 


0.6(0.7) 


0.22(0.26) 


3.3(3.9) 


Spider 


6(6) 


3.6(4.5) 


25(29) 


EBEX 


2.4(3.9) 


0.9(1.6) 


7(13) 


QUIJOTE 


0.5(1.5) 


0.34(2.0) 


2.4(6) 



TABLE III: Expected values of aBo/lO"'' Gauss for the dif- 
ferent experiments based on the experimental characteristics 
from Table II. The numbers in parenthesis corresponds to 
the results when including the foreground model, while re- 
sults from the foreground free analysis are quoted outside the 
parentheses. 



Let us first concentrate on the results obtained with- 
out the foreground model. We can identify which power 
spectrum gives the best results depends on the experi- 
ment under consideration. That TB is less important for 
QUIET than for e.g. Planck is easily understood by the 
fact that QUIET has much better Wp^ relative to w^^ 
than the bolometer based experiments. But we have also 
seen that the results from the EB and BB power spec- 
tra depend more strongly on Wp^ and w^^ than does 
the TB power spectrum. This explains e.g. why we for 
Planck find the best limits from using TB, while the best 
limits come from EB for Clover. The most impressive 
results come from the QUIET experiment's phase 2 tele- 
scopes. Here the EB power spectrum can constrain Bq 
by aso = 7 X 10~" and aso = 4 x 10"" from the 2m 
and 7m telescopes, respectively. This is approximately a 
factor 50 better than what we find for Planck TB and 
Spider, a factor 20 better than EBEX and about a factor 
5 better than Clover and QUIJOTE. 

Adding the foreground model, this picture changes sig- 
nificantly. The results from QUIET P2 now weakens by 
more than an order of magnitude, while the changes for 
the other experiments arc typically by ^ 10%. The rea- 
son why QUIET is much more sensitive to this foreground 
is caused by several factors. Firstly, the QUIET exper- 
iment covers a small fraction of the sky, and we have 
seen that while this can be beneficial in a foreground- 
free scenario, the opposite is the case when the fore- 
ground is present. Also, QUIET operates at a rela- 
tive low frequency compared to for example Clover and 
EBEX, which makes the synchrotron foreground much 
stronger, as its amplitude has a dependence. We 

have also seen that the effect of having good experi- 
mental sensitivity is less important in the presence of 
foregrounds. The result is that when including the fore- 
ground model, we find the best results from Clover EB, 
giving = 2.6 x 10~^°Gauss, which is a factor 3 bet- 
ter than what we find for QUIET, a factor 6 better than 
EBEX and an order of magnitude better than Planck. 
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Which of the results with or without the foreground 
are more close to the real performance of the experi- 
ments depends on how well the foreground can be sub- 
tracted from the maps. This will be of special importance 
for QUIET and less important for the other experiments 
considered here. Also, we have assumed perfect knowl- 
edge of the statistical properties of the foreground model. 
When analyzing a real experiment, these uncertainties 
must be taken into account and this may alter the re- 
sults, especially for a "foreground sensitive" experiment 
like QUIET. We also note that the results in Table III are 
quoted for each individual power spectrum. Simply com- 
bining the results from the different power spectra (and 
experiments) will of course improve the limits slightly. 

V. CONCLUSIONS 

In this work we have calculated CMB power spectra 
for a universe with a primordial, homogeneous magnetic 
field. We introduced a simple foreground model, and ap- 
plied a Fisher matrix approach to explore how the exper- 
imental sensitivity to such a magnetic field depended on 
different experimental parameters. We then attempted 
to forecast our ability to constrain for various future 
CMB polarization experiments. 

In our foreground model we have assumed perfect 
knowledge of the statistical behavior of the foregrounds, 
but that none of it has been removed from the CMB 
maps prior to the parameter analysis. We argue that 



this can be regarded as a conservative approach. Includ- 
ing this foreground model we find that among the ex- 
periments considered here, Clover has the best prospects 
of constraining Bq, and we find a^o « 3 x 10~^°Gauss 
in this case. When assuming that all foregrounds are 
perfectly subtracted from the maps, we find QUIET 
to be the most promising experiment, giving a^o ~ 
4 X 10~^^Gauss. This is two orders of magnitude bet- 
ter than the present upper limits on Bq from using the 
homogeneous anisotropy of CMB [9, 44]. For the QUIET 
experiment, the results rely heavily on the ability to suc- 
cessfully subtract the foregrounds, while the other exper- 
iments under consideration are less dependent on this. 

To conclude, our study shows that we will be able to 
significantly improve our knowledge of possible primor- 
dial homogenous magnetic fields in the coming years. A 
positive detection will impact both the construction of 
inflationary models and our understanding of the origins 
of magnetic fields in galaxy clusters. 
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